Walks on weighted networks 



An-Cai Wu, 1 Xin-Jian Xu, 2 Zhi-Xi Wu, 1 and Ying-Hai Wang 1 -* 
1 Institute of Theoretical Physics, Lanzhou University, Lanzhou Gansu 730000, China 
2 Department of Electronic Engineering, City University of Hong Kong, Kowloon, Hong Kong, China 

(Dated: February 6, 2008) 

We investigate the dynamics of random walks on weighted networks. Assuming that the edge's 
weight and the node's strength are used as local information by a random walker, we study two kinds 
of walks, weight-dependent walk and strength-dependent walk. Exact expressions for stationary 
distribution and average return time are derived and confirmed by computer simulations. We 
calculate the distribution of average return time and the mean-square displacement for two walks 
on the BBV networks, and find that a weight-dependent walker can arrive at a new territory more 
easily than a strength-dependent one. 
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There has been a long history of studying random 
walks to model various dynamics in physical, biological, 
social, and economic systems pj. A large body of theo- 
retical results are available for random walks performed 
on regular lattices and on the Cayley (or regular) trees 
0, , which are defined to be trees with homogeneous 
vertex degree. However, it has been suggested recently 
that more complex networks as opposed to regular graphs 
and conventional random graphs |4j are concerned to real 
worlds. Particularly, important classes of random graphs 
such as small- world networks Q and scale-free networks 
were proposed and have been examined in the last sev- 
eral years. These networks share some important prop- 
erties with real networks, such as the clustering prop- 
erty, short average path length, and the power-law of the 
vertex degree distributions. They have been applied to 
the analysis of various social, engineering, and biological 
networks including epidemic spreading, percolation, and 
synchronization, et al [71 IM la. Il0| . 

Recently, there have been several studies of random 
walks on small- world networks 0, 0, fl3L fl4| and on 
scale- free networks 0,0,0,^^. Most of studies of ran- 
dom walks focus on unweighted networks, however, the 
study of the dynamics of random walks on weighted net- 
works is missing while most of real networks are weighted 
characterized by capacities or strengths instead of a bi- 
nary state (present or absent) [HEa EJ EHH H| . In 
the weighted networks, a weight wij is assigned to the 
edge connecting the vertices i and j, and the strength of 
the vertex i can be defined as 

where the sum runs over the set v(i) of neighbors of i. 
The strength of a node integrates the information about 
its connectivity and the weights of its links. 

In this paper we study the dynamics of random walk 
processes on weighted networks by means of BBV model 



|22|. The model starts from an initial number of com- 
pletely connected vertices mo with a same assigned 
weight wo to each link. At each subsequent time step, ad- 
dition of a new vertex n with m edges and corresponding 
modification in weights are implemented by the following 
two rules: (i) The new vertex n is attached at random to 
a previously existing vertex i with the probability that is 
proportional to the strength of node i, Si/^2j Sj, imply- 
ing new vertices connect more likely to vertices handling 
larger weights, (ii) The additional induced increase 5 
in strength Si of the ith vertex is distributed among its 
nearest neighbors j 6 V(i) according to the rule 

U) ■ ■ 

viij -> + <5— , (2) 

Si 

which considers that the establishment of a new edge 
of weight wo with the vertex i induces a total increase 
of traffic S that is proportionally distributed among 
the edges departing from the vertex according to their 
weights. The BBV model suggests two ingredients of self- 
organization of weighted networks, strength preferential 
attachment and weight evolving dynamics |22j. Consid- 
ering an arbitrary finite weighted network which consists 
of nodes i = 1, . . . , N and links connecting them. The 
connectivity is represented by the adjacency matrix A 
whose element <2y = 1 if there is a link from i to j, and 
fly = otherwise. The information of edges weight is 
represented by matrix W whose element w^ is the weight 
of the edge between i and j. We restrict ourselves to 
an undirected network aij — aji and symmetrical edge's 
weight w^ = Wji. 

Assuming that edge's weight and node's strength are 
used as local information by a random walker, we define 
two kinds of walks, weight-dependent walk and strength- 
dependent walk. For weight-dependent walk, a walker 
chooses one of its nearest neighbors with the probability 
that is proportional to the weight of edge linked them. 
The transition probability from node i to its neighbor j 
is 
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When time becomes infinite, one can find the walker stay- 
ing at node i with the probability Pf, which is defined 
as the stationary distribution. Following the ideas devel- 
oped by Noh and Rieger [l^, we can write 



pu 



(4) 



namely, the larger strength a node has, the more often it 
will be visited by a random walker. 

For strength-dependent walk, a walker at node i at 
time t selects one of its neighbors with the probability 
which is proportional to the selected node's strength to 
which it hops at time t + 1. The transition probability 
from node i to its neighbor j is 



Pf 



Sj dij 



(5) 



where s i — Y^i^vli) Sl an< ^ v tf) denotes the set of all 
neighboring vertices of node i. Supposing that a walker 
starts at node i, and the probability that the walker at 
node k after t time steps denoted by P^ k {t), then the mas- 
ter equation for the probability -P/_>j to find the walker 
at node j at time t + 1 is 

/(t+1) = £-a* 



p? 



(6) 



An explicit expression for the transition probability 
P*_>j(t) to go from node i to node j in t steps follows 
by iterating Eq. JfjJ 
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(7) 



Comparing the expressions for P-^Jt) and P^^t), we 
get 



SiSiPl^At) 



s i s jPj- 
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(8) 



This is a direct consequence of the undirectedness of the 
network. We can also define the probability Pf as the 
stationary distribution when the evolving time becomes 
infinite. Eq. (8J) implies that SiS^Pj — s.,s„P/, and there- 
fore one can obtain 0] 



'j^- 1 i > 



SiS, 



(9) 



Now we discuss the stationary distribution for edges, 
i.e., the probability that an edge is chosen by the walker 
to follow as the evolving time becomes infinite. In un- 
weighted networks, all edges are equal and a random 
walker will choose one of its neighboring edges at the 
same probability. So each edge in the network has the 
same probability to be chosen by the walker when t — > oo 
for weight-dependent walk. In weighted networks, how- 
ever, the walker will choose an edge according to the 
weight of it or the strength of the node connected by 
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FIG. 1: Log-log plots of P™ vs. s< (a) and P e ™ vs. tu y (b) in 
the BBV network with N = 1000, mo = 3 and different values 
of S. Slopes of all the curves are equal to 0.9999 ±0.0006. The 
data were obtained after walking 10 8 steps on the network. 



it. Then the relation between the stationary distribution 
for edges P eij and the stationary distribution for nodes 
Pi can be written as 



P^.aa Pi Pi- 



(10) 



where Pi^j is the transition probability from node i to 
nodej. For weight-dependent walk, substituting Eqs. © 
and into Eq. (|10|) , we obtain the stationary distribu- 
tion for edges 



2wi 
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(11) 



For strength-dependent walk, substituting Eqs. © and 
© into Eq. ljTU|) . we obtain the stationary distribution 
for edges 



ps 



(12) 



In Fig. Owe plot P t w vs. (Fig. Ufa)) and P e ™ vs. 
Wij (Fig. IHb)) in log- log scale in the BBV network. The 
power-law property of Eqs. © and l|ll|l is presented in 
excellent agreement with the numerical results. In Fig. 
|2]we show the log-log plots of P* vs. SjS, (Fig. Ufa)) 
and P e s .. vs. s.Sj (Fig. ©b)) in the BBV model. The 
numerical results are also in good agreement with Eqs. 
© and (El- 
Next we study average return time which is the average 
time spent by a walker to return to its origin. From its 
definition, we can easily obtain that the average return 
time is equal to the reciprocal of the stationary distribu- 
tion. The average return time for node i is 



Ej «j 



(13) 
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FIG. 2: Log-log plots of P" vs. SiSi (a) and P".. vs. SiSj (b) in 
the BBV network with N = 1000, mo = 3 and different values 
of S. Slopes of all the curves are equal to 1.018 ± 0.003. The 
data were obtained after walking 10 s steps on the network. 



for weight-dependent walk and 
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for strength-dependent walk, respectively. 

Using the same methods, the average return time for 
edge ey can also be obtained 



CO 
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for weight-dependent walk and 
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(15) 
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for strength-dependent walk, respectively. 

In Fig. |31 we show the log-log plots of (T$) vs. s(i) 
and (TS) vs. s(i) in the BBV networks. The slope of 
the curve in Fig. 02a) is -0.9903 ± 0.0007, consistent 
(|13H . The BBV networks have two proper- 
(i) node's strength is proportional to node's 
~ ki\ (ii) the average nearest neighbor de- 
(fc) ~ fc-2+V/3 with p = (25 + i)/(2J + 2). 



with Eq. 
ties in [2i 
degree, s 
gree is k, 



Considering these two ingredients, one can observe that 
s i s. l ~ s 1 / 13 , and obtain (TS) ~ s i l ^ . Fig. [Jb) gives 
that — 1//3 — —1.355 ± 0.006 which agrees with the the- 
oretical value -1/(3 = -{25 + 2)/ (25 + 1) = -4/3. In 
Fig. |3 the slope of the curve for strength-dependent 
walk is steeper than that for weight-dependent walk, giv- 
ing rise to a broader distribution of average return time 
for strength-dependent walk. In order to confirm this 
point, we derive the expression for the distribution of 
average return time. In BBV network, we know that 
the strength distribution behaves as P(S) ~ S* -7 , where 
7 = (3 + 45)/ (1 + 25) H3- According to the above analy- 
sis, we can obtain that the distribution of average return 



FIG. 3: Numerical results for the average return time of 
node for the two different walks , weight-dependent walk (a) 
and strength-dependent walk (b), on the BBV network with 
N — 1000, mo = 3 and 5 = 1.0. The data were averaged on 
1000 networks and obtained after walking 10 6 steps in each 
network. 



time is P{T W ) ~ (T u, )- 1 /(i+<5) f or t h e weight-dependent 
walk, and P(T S ) ~ 1.0 for the strength-dependent case. 
Thus, we can see, for those nodes with large strength 
the strength-dependent walker spends more time in vis- 
iting them than that the weight-dependent walker does. 
This point can be reflected by the difference of the mean- 
square displacement for two walks which is shown in the 
following. 

The mean-square displacement (R 2 ) is a measure of 
the distance R covere d by a typical random walker after 
performing t steps [13, E3 • To calculate this quantity we 
first, at each time step, find the minimal distance from 
the current position of the walker to the origin (i.e., the 
smallest number of steps needed for the walker to reach 
the origin ) using a breadth-first search method. Then 
we allow the walker to move through the network until 
< R 2 > has saturated. Finally, we average over different 
initial positions of the walker and realizations of the net- 
work. Fig. 0] shows the results of (R 2 ) as a function of 
time t. We note that (R 2 ) equilibrates after a few steps 
to a constant displacement value. This is a simple mani- 
festation of the very small diameter of the BBV network. 
Note also that the value of plateau is higher for weight- 
dependent walk than that for strength-dependent walk. 
That is, a strength-dependent walker spends more time 
in visiting these nodes with large strength or large de- 
gree, i.e. a weight-dependent walker can arrive at a new 
territory more easily than a strength-dependent walker. 

In summary, we have studied two different types of 
walks on weighted networks, weight-dependent walk and 
strength-dependent walk. We derived exact expressions 
for the stationary distribution and the average return 
time for the two walk processes, and confirmed them by 
simulations on BBV networks. Then we analysed the dis- 
tribution of average time for two walks and found that 
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it is broader for dependent-strength walk than that for 
dependent-weight walk on the BBV network. Finally we 
computed the mean-square displacement (R 2 ). For both 
walks, (R 2 ) was found to reach the saturation after a 
few time steps which is a result of the very small di- 
ameter of the underlying graph. Furthermore, the differ- 
ence of average-square displacement for two walks implies 
that a weight-dependent walker can arrive at a new ter- 
ritory more easily than a strength-dependent walker on 
the BBV network. 



FIG. 4: Mean-square displacement (R 2 ) as a function of time 
t for two kinds of walks, weight-dependent walk (closed cir- 
cles) and strength-dependent walk (closed squares), on the 
BBV network with N = 10 5 , m = 3, and 5 = 1.0. All the 
plots were averaged over 1000 different initial positions of the 
walker and 20 realizations of the network. 
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